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Abstract: We present the differential rates and branching ratios of the radiative decays 
r —)■ with I = e or pL, and /i —?■ in the Standard Model at next-to-leading order. 

Radiative corrections are computed taking into account the full depencence on the mass 
mi of the final charged leptons, which is necessary for the correct determination of the 
branching ratios. Only partial agreement is found with previous calculations performed 
in the m; —)• 0 limit. Our results agree with the measurements of the branching ratios 
B{fi —)• euuy) and B{t —?■ pii'iz'y) for a minimum photon energy of 10 MeV in the pL and 
r rest frames, respectively. Babar’s recent precise measurement of the branching ratio 
B{t ^ eHiz'y), for the same photon energy threshold, differs from our prediction by 3.5 
standard deviations. 
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1 Introduction 

Muon and r leptonic decays offer one the most powerful tools to study the Lorentz structure 
of weak interactions. Their theoretical formulation via the Bouchiat-Michel-Kinoshita- 
Sirlin (BMKS) parameters [1-4] places them in a unique position to investigate possible 
contributions beyond the V-A coupling of the Standard Model (SM). Radiative /i and r 
leptonic decays, where an inner bremsstrahlung photon is measured, can be predicted with 
very high precision and provide an independent determination of the BMKS parameters 
as well as the possibility to extract new combinations like the fj parameter [5-7]. A new 
preliminary measurement of the muon ff was reported recently [8] , while analyses are ongoing 
to determine the 77 and parameters of the r [9, 10]. Precise data on radiative r leptonic 
decays also offer the opportunity to probe the electromagnetic properties of the r and 
may allow to determine its anomalous magnetic moment which, in spite of its precise SM 
prediction [11], has never been measured [12-14]. 

Recently, the Babar collaboration performed the measurements of the r —>■ (I = 
e, 77 ) branching fractions for a minimum photon energy cjq = 10 MeV in the r rest frame [15, 
16]. The experimental precision of these measurements, around 3%, requires the SM predic¬ 
tion of the branching ratios at next-to-leading order (NLO). Indeed these radiative correc¬ 
tions are not protected from mass singularities by the Kinoshita-Lee-Nauenberg (KLN) the¬ 
orem [17-19] and are therefore expected to be of relative order (a/vr) ln(mi/mr) ln(wo/mr), 
corresponding to a large 10% correction for I = e, and 3% for I = /i. Furthermore, spe¬ 
cial attention must be paid to the role played by the final lepton mass mi and the limit 
mi —7- 0. In fact, in apparent contradiction to the naive expectation based on the V-A weak 
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interaction and helicity conservation in massless QED, final state charged leptons have a 
finite probability of being right-handed even in the chiral limit m; —?■ 0. This non-intuitive 
feature is a consequence of helicity-flip bremsstrahlung in QED and appears as a peculiar 
mass-singularity cancellation in the collinear region [19-22], 

Radiative /U and r leptonic decays also constitute an important source of background 
for experiments searching for charged lepton flavour violating decays, such as —)• 

T —>■ Ij, and even /r'*' —>■ because of the internal conversion of photons to electron- 

positron pairs. In the next stage of the Meg experiment [23], as well as in future searches 
for —)> e'’'e“e'’' at the Mu3e experiment [24], the desired sensitivity will require a refined 
control and a precise measurement of these backgrounds, which are indistinguishable from 
the signal except for the missing energy carried away by the neutrinos. An improvement of 
their theoretical calculation down to a precision of 0(1%) is desired [25]. Preliminary new 
measurements of the branching ratio [i —)• were presented recently by the Meg [25] 

and PiBETA [8] collaborations. 

We begin our analysis in section 2 reviewing the SM prediction for the differential 
decay rates at leading order (LO). The decay rates at NLO are presented in section 3. 
Our NLO predictions for the branching ratios of radiative /r and r leptonic decays are 
reported in section 4, where they are compared with published experimental measurements. 
Conclusions are drawn in section 5. 


2 Differential decay rates at LO 

The LO SM prediction for the differential rates of the radiative decays 

( 2 . 1 ) 
( 2 . 2 ) 


/i e 

T~ r n Vr 1, 


with / = e or /i, of a polarized ^ or r in their rest frame is 

d^PLo aGlM^ xl3 

_ — _ £ _ 

dx dy dQi dQ-y (47r)® l + dwim^jUie) 


GLoix,y,c) + xl3h-piJi^o{x,y,c) + yn ■ Py Ki^oix,y,c) 


, (2.3) 


where Gp = 1.166 378 7(6) x 10“^ GeV“^ [26] is the Fermi constant, defined from the muon 
lifetime, and a = 1/137.035 999157 (33) is the fine-structure constant [27, 28]. Calling M 
and m the masses of the initial and final charged leptons (neutrinos and antineutrinos are 
considered massless) we define r = m/M and rw = Af/Mw, where Mw is the VP-boson 
mass; p and n = (0, h) are the four-momentum and polarization vector of the initial p or r, 
with = —1 and n-p = 0. Also, x = 2£';/M, y = 2ujIM and f3 = \pi\/Ei = y^l — 
where pi = (Eppi) and py = {uj,py) are the four-momenta of the final charged lepton and 
photon, respectively. The final charged lepton and photon are emitted at solid angles D/ 
and D,),, with normalized three-momenta pi and py^ and c = cos0 is the cosine of the angle 
between pi and py. The corresponding formula for the radiative decay of a polarized or 
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//■*■ is simply obtained inverting the signs in front of the scalar products h • pi and h ■ in 
eq. (2.3). The term 6w{m^,me) = 1-04 x 10“® is the tree-level correction to muon decay 
induced by the IT-boson propagator [29, 30]. 

The function Glo and, analogously, Jlo and Ki^q, is given by 


Gi^o{x,y,c) 


4 


go{x,y,z) + r‘^g^{x,y,z) + O (r^) , 


(2.4) 


where 2 ; = xy {1 — c/3) /2. The functions go, jo, and ko, computed in refs. [31-35], arise 
from the pure Fermi V-A interaction, whereas g^/, j^, and fcwi calculated in ref. [30], 
are the leading contributions induced by the VF-boson propagator. For the decay (2.1), 
~ 2 X 10“®, while for (2.2), ~ 5 x 10 Distributions in terms of the helicities of the 

final lepton and photon were studied in ref. [36]. If the initial or are not polarized, 
eq. (2.3) simplifies to 


d^FLo aG'^pM^ Svr^ x/3 

dxdcdy (dvr)® l + 6w{m^,me) 


3 Differential decay rates at NLO 

We will now consider the SM prediction for the differential rates of the decays (2.1,2.2) at 
NLO in a. These NLO corrections were computed using the effective Fermi Lagrangian, 
i.e. collapsing the SM weak decay vertices, mediated by the LF-boson, to a four-fermion 
interaction. In this approximation, tiny terms of 0{arnj^/M^) ~ 10“® and 0{am‘^/M^) ~ 
10“® were neglected, but they are expected to be even smaller than the uncomputed next- 
to-next-to-leading order (NNLO) corrections of Throughout the calculation, the 

full dependence on the mass ratio r = m/M has been taken into account. 

3.1 QED radiative corrections 

In the Fermi V~A theory, a virtual photon can only be exchanged between charged fermions; 
see figure 1. We performed the computation of the one-loop QED diagrams via the standard 
Passarino-Veltman reduction of the tensor integrals [37] using FORM [38] and the Mathemat- 
ica package FeynCalc [39]. All IR-convergent scalar integrals were computed analytically 
following [40] and checked numerically with LoopTools [41]; the IR divergent ones were 
taken from refs. [42, 43]. 

UV-finite results were obtained in the on-shell renormalization scheme. Indeed, as 
shown long ago by Berman and Sirlin [44], to leading order in Gjr, but to all orders in a, the 
radiative corrections to muon decay in the Fermi V~A theory are finite after mass and charge 
renormalization. A small photon mass A was introduced to regularize the IR divergences, 
while the mass ratio r regularized the collinear ones. As a check of the calculation of the 
renormalized one-loop amplitude e*^{p-y)Ai^i^^, where e*^{p^) is the polarization vector of the 
outgoing photon, we explicitly verified that satisfies the Ward identity [p-y)^ = 0- 


Experimentally, double bremsstrahlung events in which one of the two photons is too 
soft to be detected are counted as single-photon radiative decays (2. 1,2. 2). The "soft" 
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Figure 1. Radiative and r leptonic decays: one-loop QED corrections. 


differential decay rate for these events should therefore be added to the LO differential 
decay rate in eq. (2.3) and to the virtual QED corrections discussed above, thus removing 
the IR divergence of the latter. 

Specifically, let us consider decay events with the emission of two photons with energies 
Lj and uj' (and normalized energies y = 2u}jM and y' = 2ui'/M), and let us assume that 
photons can be measured when their normalized energies are above a given threshold y^. 
In the limit <C 1, the differential decay rate for the double photon emission obtained 
integrating y' up to y^ is [17, 45] 


(y'; 


dx dy dQ .1 dVl^ 


= -r-| (Ms/D" - In 4, 


+ iLi'2 


vr 

2/3 

1 + /3 


--r In 

2/3 


M2 

1 + /3 




1 + /3 
1-^ 
1 + /3 

1-/3 


+ 


- 1 


d^T 


LO 


dx dy dO^i dil^ 


(3.1) 


with y > yQ. We checked that the IR divergence of this soft bremsstrahlung contribution is 
canceled by that arising from the virtual QED corrections. 

In order to calculate the branching fractions of the radiative decays (2.1,2.2) at NLO 
(see section 4) , we also computed the differential rate of the leptonic decay of an unpolarized 
r or /i with "hard" double-photon emission, which was obtained integrating y' from the 
threshold y'^ up to its kinematic upper limit. The calculation of this hard bremsstrahlung 
differential rate required a detailed analysis of the allowed five-particle phase space which 
was then integrated numerically with the CUBA library [46]. 


3.2 NLO results 

The NLO prediction for the differential rate of the radiative fj, and r leptonic decays 
(2. 1,2. 2), with the possible emission of an additional soft photon with normalized energy 
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y' < y'o, is 

d^T{y'o) aGlM^ 


X 


dx dy dQi dQ-y (47r)® l + 6w{'m^,me) 


X G + xP h -Pi J + yh-pyK + xy/d h ■ {pi x py) L , (3.2) 


with y > yg. The function G(x,y,c;yg) and, analogously, J and K, is given by 



where gQ{x,y,z) and g^fir{x,y,z) are the LO contributions described in section 2, while 


fl'NLo(3^, y, 2 ; yg) is the sum of the virtual and soft bremsstrahlung contributions (3.1). 


The function L{x,y,z), appearing in front of the term h ■ {pi x py), is only induced by 
the loop corrections and is therefore of O^a/n). In particular, L{x,y,z) is of the form 
Pnix, y, z) Im [In{x, y, z)]^ where Pn are polynomials in x, y, z and Inix, y, z) are scalar 
one-loop integrals whose imaginary parts are different from zero. The functions G, J, K 
and L are free of UV and IR divergences. Their (lengthy) explicit expressions are provided 
in attachment to this paper. ^ If the initial or are not polarized, eq. (3.2) simplihes 


to 


d^T (yg) aG‘^pM^ Svr^ xjd 

dxdcdy (dvr)® 1 + 


Gix,y,c] y'o). 


(3.4) 


QED one-loop corrections to radiative muon decays were previously computed in refs. [34, 
45] at different levels of completeness, or as part of NNLO corrections to muon decay [47- 


49]. The isotropic correction yNco was computed by the authors of ref. [34], but we couldn’t 
get hold of the hies with their expressions. The authors of ref. [45] performed the calcu¬ 
lation of radiative p decays with the full spin dependence, but in the r = rue/m^ —)• 0 
limit, ^ whereas, as we will discuss in section 4, terms in G proportional to cannot be 


neglected if the differential decay rates are integrated over the entire allowed phase space. 
We compared our isotropic function yNcO) taking for this test the limit r —?■ 0, with the 
result yNLo|r-!>0 of ref. [45], hnding perfect numerical agreement. On the contrary, our re¬ 
sults for the anisotropic contributions Jnlo and fc^Lo differ, even in the r —)■ 0 limit, from 
those of ref. [45]. Moreover, the contribution of the function L{x,y,z) has been previously 
overlooked. 

4 Branching Ratios 

The kinematic limits of integration for the variables x, c, and y are 


2r<x<l + r‘^, -l<c<l, 0 < y < ymax(a;, c), 


(4.1) 


where the maximum normalized photon energy as a function of x and c is 



(4.2) 


Fortran code with all the LO and NLO contributions discussed in this paper is available upon request. 
^This limit is appropriate if the collinear region is excluded from the phase space Integration. 
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However, every experimental setup has a minimum photon energy wq = yo{M/2) below 
which photons are not detected. As the constraint uq < ymax(a^;C), necessary to measure 
radiative decays, leads to the bound c < Cmax(a^)) with 


Cmax(3^) — 


2 (l + - x) - (2 - x)yo 

X13 yo 


the kinematic ranges of x, c, and y > yo are reduced to 


(4.3) 


2r < X < 1 + r^, -1 < c < min{l, Cmax(3;)}, yo < y < ymax(x, c). (4.4) 


Integrating the differential rates in eq. (3.4) over the kinematic ranges (4.4) and multiplying 
them by the ^ or r lifetimes = 2.1969811(22) x 10“® s and Tr = 2.903(5) x 10“^^ s [26], 
we obtain the NLO predictions for the branching ratios of the radiative decays (2.1,2.2). 


4.1 Branching ratios at LO 

The analytic integration over the kinematic ranges (4.4) of the LO differential rate in 
eq. (2.5) with a minimum energy yo = 2ojo/M gives [31, 33] 

Llo 




^ im) = 3Li2(yo) - y + (inr + L) (61nyo + 6yo + yo) + 

+ ^ (125 + 452/0 - 332/0 + ^vl) 2 /o - ^ (6 + 2 /o) yo Inyo, 


(4.5) 


(4.6) 


where yo = l~yo nnd the dilogarithm is defined by Li 2 (z) = — j^dt Terms depending 

on the mass ratio r have been neglected in the expression for Lf(yo), with the obvious 
exception of the logarithmic contribution which diverges in the limit r —)■ 0. However, 
as already mentioned in section 3, terms in the integrand Gi,o{x,y,c) proportional to 
were not neglected when performing the integral to obtain (4.6), as they lead to terms of 
0(1) in the integrated result H(yo). In fact, the functions multiplying these terms in 
the integrand generate a singular behavior in the r —)• 0 limit after the integration over 
c = cosO: terms proportional to j in Gi^o{x,y,c) lead to a nonvanishing contribution 
to the integrated decay rate since fdc(l/z^) oc 1/z \s evaluated at the integration limit 
c —>■ 1 where z —)> xy (1 — /3)/2 ~ r‘^{y/x) for x 3> 2r. This feature, first noted in [19], is due 
to the appearance of right-handed electrons and muons in the final states of (2.1,2.2) even 
in the limit r —>■ 0, and is a consequence of helicity-flip bremsstrahlung in QED [19-22]. We 
also note that the presence of the mass singularity Inr in the integrated decay rate Tlq (yo) 
does not contradict the KLN theorem, which applies only to total decay rates [17-19]. The 
tiny corrections induced by the LL-boson propagator were neglected in eq. (4.5). 

If we multiply the analytic result for Tlq (yo) in 6<1- (4-5) by the lifetimes with a 
threshold cjq = 10 MeV we obtain the following LO predictions for the branching ratios: 
1.83 X 10“^ (r —> eT'i^y), 3.58 x 10“^ (r —^ yi^z^7), and 1.31 x 10“^ (// —> eiyuj). These 
values are in perfect agreement with the results Hlo reported in table 1, obtained integrating 
numerically the LO differential rates in eq. (2.5), with the exception of the r —>• /rPz^y 
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value; this difference is due to the terms neglected in the analytic result (4.5). We note 
that, had we performed the numerical integrations of the LO differential rates neglecting 
terms proportional to in the integrand Gi^o{x, y, c), we would have obtained the incorrect 
results 1.96 x 10“^ (r —)■ er'z^y), 4.92 x 10“^ (r —)■ and 1.44 x 10“^ (/r —)■ euv'^). 


4.2 Branching ratios at NLO 

The branching ratios of the radiative decays ( 2 . 1 , 2 . 2 ) can be distinguished in two types: 
the "exclusive" ones, {vo): ^re measured counting the decays with the emission of one, 
and only one, photon of energy larger than the detection threshold yg- Such measurements 
clearly include double bremsstrahhmg events in which one of the two photons is too soft to 
be detected, but decay events with the emission of two hard photons, both with energies 
larger than yo; must be discarded. The "inclusive" branching ratios (yg) are dehned 
measuring radiative decays with the emission of at least one photon with energy higher 
than yg. In these inclusive measurements, if two photons with energies y and y' larger than 
yg are detected, the event must be included. 

At LO, the theoretical predictions for these exclusive and inclusive branching ratios 
clearly coincide - double bremsstrahlung events are simply not considered. At NLO, the 
theoretical prediction for the exclusive branching ratios can be calculated integrating the 
differential rate d?T (yg) in eq. (3.4), multiplied by the lifetime, setting yg = yg. As d^T (yg) 
was obtained adding the analytic formulae for the LO contribution, the virtual corrections, 
and the soft term (yg) (see eq. (3.1)), thus cancelling the IR divergences, this result 

will be called 

^Lvs(yo,yo) = ^Lo(yo) + ^vs(yo, y'o), (4-7) 

with yg = yg, where the LO term has been separated for convenience.^ However, as 
d^T^f (y') was calculated analytically in the limit yg <C 1, Hlvs(2/0) 2/0 ) is inadequate to pre¬ 
dict the exclusive branching ratios if yg is large. This is the case, for example, for y —)• 
decays with cjg = 10 MeV (yg ~ 0 . 2 ), where the NLO correction Bvs{yo,yo) = — 1-1 x 10 “^ 
obtained in this manner turns out to be off by a factor of two w.r.t. the correct value 
— 2.2 X 10 “^ (see ( 2 / 0 ) later on). We will therefore proceed differently, deriving the 
exclusive branching ratios from the inclusive ones. 

The inclusive branching ratios can be calculated adding to Hlvs( 2/0) 2/g) the branching 
ratios for the hard double bremsstrahlung events. More precisely, let us call Hi^^’^'^(yg, yg) 
the branching ratio of hard double bremsstrahlung decays where both photons have energies 
larger than yg, and at least one of them has energy above yg. The NLO prediction for the 
inclusive branching ratio is given by 


S'""(2/o) 


lim 

y'o^o 


Slvs(2/o,2/o) + ^ 


hard 

77 


( 2 / 0 , 2 / 0 ) 


(4.8) 


The sum in square brackets of eq. (4.8) should not depend on yg, because the sum of soft 
and hard double-photon emission describes events with one photon of energy larger than yg 
plus another photon of any energy allowed by the kinematics. Therefore, the yg-dependent 


®Note that, in general, the two thresholds yo and j/g can be different, with y), < yo- 





terms of the soft radiation contribution should be canceled by those of the hard one. This, 
however, is not exactly the case because, as mentioned earlier, the soft term was calculated 
in the <C 1 limit and the cancellation of the j/g-dependent terms is therefore not exact. 
For this reason we introduced the limit j/g —?• 0 in eq. (4.8). 

Once the inclusive branching ratios are calculated, the exclusive ones can be computed 
subtracting 

13^- (yo) = (yo) - (yo, yo) • (4.9) 

The second term on the r.h.s. of eq. (4.9) subtracts in fact those events in which both photons 
have energies above the threshold yg. Such events are discarded in exclusive measurements 
of the branching fractions. We point out that B^^‘^{yo) differs, in general, from the branching 
ratio ;Blvs( 2 / 0 ) 2/o) discussed above (we already mentioned the discrepancy for radiative y 
decays with coq = 10 MeV). For convenience, B^^'^{yo) and .B^'^'^(yo) are separated into LO 
contributions and NLO corrections: 

B^^^ (yo) = B^o {yo) + B^ro (yo) , (4.10) 

(yo) = B^o (yo) + BiTo (yo) • (4.ii) 


4.3 Numerical results 

Exclusive and inclusive branching ratios for the radiative decays (2.1,2.2) are presented 
in table 1 for a threshold wq = yo (-4^/2) = 10 MeV. All branching ratios were computed 
keeping into account the full dependence on the mass ratio r. In fact, also at NLO, for 
the same reasons discussed in the LO case, terms in the integrand G(x, y, c) proportional 
to cannot be neglected when performing the integrals to calculate the branching ratios. 
Numerical integrations were performed with the Cuba library [46] and tested with differ¬ 
ent numerical integration methods. The hard double bremsstrahlung terms .B]^“^(yo, yg), 
necessary both for .B^™(yo) and for B^°^(yo), required a detailed study of the five-particle 
phase space which was then integrated numerically. The values obtained for y^)) 

were checked with MadGraphS [50]. 

Uncertainties were estimated for uncomputed NNLO corrections, numerical errors, and 
the experimental errors of the lifetimes. The former were estimated to be ~ 

(a/vr) Inr ln(a;o/M) For cvg = 10 MeV they are about 10%, 3%, and 3% for 

r — )• T —7- and y — )• respectively. They appear with the subscript 

"N" in table 1. The branching ratios due to hard triple-photon emission, estimated with 
MadGraph5, are much smaller than these uncertainties. Numerical errors, labeled in ta¬ 
ble 1 by the subscript "re", are smaller than those induced by missing radiative corrections. 
These two kinds of uncertainties were combined to provide the theoretical error of the final 
and .6^“ predictions, labeled in table 1 by the subscript "th". The uncertainty due 
to the experimental error of the lifetimes is labeled by the subscript "r"; it is negligible in 
radiative y decays. 

The recent measurements by the Babar collaboration of the branching ratios of the 
radiative decays r — ?■ /rerey, with I = e and y, for a minimum photon energy ujq = 10 MeV 
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r —^ ei'U'y 

r —)• niyi"y 

—)• evi'j 

Bpo 

1.834 X 10-2 

3.663 X 10-3 

1.308 X 10-2 

r?lnc 

^NLO 

-1.06(1)„(10)7v X 10-3 

-5.8(1)„(2)^ X 10-3 

-1.91(5)„(6)7v X 10-^ 

/^Exc 

^NLO 

-1.89(1)„(19)7v X 10-3 

-9.1(1)„(3)7v X 10-3 

-2.25(5)„(7)7v X 10-4 

^Inc 

1.728 (10)th(3), X 10-2 

3.605 (2)th(6)r X 10-3 

1.289 (l)th X 10-2 

^Exc 

1.645 (19)th(3)r X 10-2 

3.572 (3)th(6)r X 10-3 

1.286 (l)th X 10-2 

Bexp 

1.847 (15)gt(52)sy X 10-2 

3.69(3)gt(10)sy X 10-3 

1.4(4) X 10-2 


Table 1. Branching ratios of radiative ^ and t leptonic decays with minimum photon energy loq = 
10 MeV. Inclusive and exclusive predictions are separated into LO contributions (Slo) 

and NLO corrections Uncertainties were estimated for uncomputed NNLO corrections 

(JV), numerical errors (n), and the experimental errors of the lifetimes (r). The first two types of 
errors were combined to provide the final theoretical uncertainty (th). The last line reports the 
experimental measurements of refs. [15, 52]. 


in the r rest frame, are [15, 16]: 

Sbxp (t ^ ei>i^7) = 1-847 (15)st(52)sy x lO'^, (4.12) 

^EXP (r ^ = 3-69 (3)st(10),y x lO’^, (4.13) 

where the first error is statistical and the second is systematic. These results are sub¬ 
stantially more precise than the previous measurements 1.75(6)gt(17)sy x 10“^ (r —)■ eT'i^y) 
and 3.61(16)st(35)sy x 10“^ (r —)• of the Cleo collaboration [51]. The signature 

for r —)• Wi"y decays is a charged particle and a photon. In particular, in the Babar 
measurements each e'’'e“ —)• event was divided into hemispheres (signal and tag) in 

the CM frame. Depending on the signal mode, either a muon or an electron with a single 
photon candidate were required on the signal side [15, 16]. The experimental values in 
eqs. (4.12,4.13) should therefore be compared with our predictions for the exclusive branch¬ 
ing ratios 1.645 (19)th(3)r x 10“^ and 3.572 (3)th(6)r x 10“^, respectively (see table 1). 
For T —)• decays, the branching ratio measurement and prediction agree within 1.1 

standard deviations (l.lu). On the contrary, the experimental and theoretical values for 
r —>■ decays differ by 2.02 (57) x 10“^, i.e. by 3.5cr. This puzzling discrepancy deserves 

further researches. From table 1 we also note that the inclusive branching ratio of r —>■ 
is less sensitive to the NLO corrections than the exclusive one. 

The branching ratio of the radiative decay /i —>■ ezzi/y was measured long ago for a 
minimum photon energy ujq = 10 MeV in the /U rest frame [52], 

-Bexp {b = 1.4 (4) X 10“^. (4-14) 

This measurement agrees with our theoretical prediction, and new precise results are ex¬ 
pected to be published in the near future by the Meg and Pibeta collaborations [8, 25]. 
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5 Conclusions 


In this work we studied the SM prediction of the differential rates and branching ratios of the 
radiative decays r —?■ (I = e, and /i —)• The NLO corrections were computed 

using the effective four-fermion Fermi Lagrangian plus QED, taking into account the full 
depencence on the mass ratio r = m/M. The resulting differential rates for the radiative 
decays of a polarized /r or r were presented in section 3, eq. (3.2). There, the functions G, 
J, K, and L contain the LO contributions, inclusive of tiny effects induced by the IT-boson 
propagator, the virtual corrections at NLO, and the "soft" double bremsstrahlung decay 
rate due to events in which one of the two photons is too soft to be detected. Explicit 
analytic expressions for G, J, K, and L are provided in attachment to this paper. The 
differential rate of the leptonic decay of an unpolarized or r with "hard" double-photon 
emission, necessary to evaluate the branching ratios, was also calculated. Agreement was 
found with an earlier calculation of the isotropic function G, which was however performed 
in the r —)> 0 limit. On the contrary, our anisotropic J and K functions differ from an 
earlier calculation, even in the massless r —)> 0 limit, and the anisotropic function L has 
been previously overlooked. 

The branching ratios for a minimum photon energy cjq = yo{M/2) were presented in 
section 4 integrating the differential rates of section 3 over the allowed kinematic ranges. 
Particular attention was paid to terms proportional to in the integrand G, as they lead 
to terms of relative 0(1) in the integrated result - they cannot be neglected. Branching 
ratios were distinguished in "exclusive", (uo): measured counting the decays with the 
emission of one, and only one, photon of energy larger than the detection threshold i/q, 
and "inclusive" ones, (uo), defined measuring radiative decays with the emission of at 
least one photon with normalized energy higher than i/q. Precise numerical results were 
presented for both of them, in table 1, for ujq = 10 MeV. Uncertainties were estimated 
for uncomputed NNLO corrections, numerical errors, and the experimental errors of the 
lifetimes. 

Our predictions agree with the old measurement of the branching ratio B{y ^ 
and precise new results for this decay are expected to be published in the near future by 
the Meg and Pibeta collaborations. Also the recent precise measurement by Babar of 
the branching ratio B(t —)• /rPz/y), for too = 10 MeV, agrees with our prediction within 1.1 
standard deviations (l.lcr). On the contrary, Babar’s recent measurement of the branching 
ratio B(t —>■ ePi^y), for the same threshold wq, differs from our prediction by 3.5a. This 
puzzling discrepancy deserves further researches. 
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